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Abstract: We consider a disordered two-dimensional electronic system in the limit of high 
magnetic field at the metal-insulator transition. Density of states close to the Fermi level 
acquires a divergent correction to the lowest order in electron-electron interaction and shows 
a new power-law dependence on the energy, with the power given by the anomalous diffusion 
exponent rj. This should be observable in the tunneling experiment with double-well GaAs 
heterostructure of the mobility ~ 10 4 cm 2 V/s, at fields of ~ 10T at temperatures of ~ 10mA" 
and voltages of ~ lfiV. 

PACS: 71.30+h, 73.40 Hm, 73.40 Gk 



1 Introduction 



One-parameter scaling theory of localization predicts that, at least in the absence of time- 
reversal breaking fields and interactions, there is a single independent critical exponent, 
the one for the localization length, which characterizes the metal-insulator transition as a 
critical point In particular, the exponent which governs the diffusion behavior and the 
eigenstates correlations at the mobility edge is set by one-parameter scaling hypothesis to 
i] = 2 — d, where d is the dimensionality of the system. In strong magnetic field however, 
the observation of integer quantum Hall effect signals a breakdown of one-parameter scaling 
in two-dimensional (2d) electron system 0. Within each disorder-broadened Landau level 
(LL) there is an unique value of the energy where the extended states reside and the system 
undergoes a metal-insulator transition when the Fermi level is at the mobility edge. Chalker 
and co-workers |§ demonstrated another surprising feature of the localization transition in 
strong fields: right at the mobility edge the diffusion constant develops a dependence on the 
frequency and the wave-vector, D(q 2 /u) oc (q 2 /uj)~ t 'I 2 for large q 2 /uj, with a non-zero value 
of the exponent rj. Equivalently, the wave-functions at the mobility edge are neither truly 
extended nor localized, but instead are fractals with the generalized dimension <i e // = 2— 77 [Q 
||. The value of the exponent 77 is expected to be an universal quantity which characterizes 
the strong-field limit and independent of the details of a disorder potential ||. 

The result that the exponent 77 is non-trivial in high magnetic field has been conjectured 
from an exact inequality satisfied by the two-particle spectral function in the lowest LL 
and by using the hypothesis of scale-invariance at the mobility edge 0. Subsequently, the 
conclusion has been confirmed in numerical calculations @ , [[J , || • Experimentally however, 
the anomalous diffusion exponent has yet to be observed, and in that context it is important 
to find measurable consequences of this new scaling behavior at the mobility edge. Shimshoni 
and Sondhi |7] have argued that the temperature dependence of Coulomb drag between two 
electron layers at the transition yields the information on the exponent 77. Brandes et al. 
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H have suggested that r\ should appear in the power of the temperature dependence of the 
energy loss rate of hot electrons due to phonon emission and of the quasiparticle lifetime due 
to electron-phonon interaction. The quasiparticle lifetime determined by electron-electron 
interactions however, is not altered by the anomalous diffusion exponent in any significant 
way |§, |§. In this paper we demonstrate that when the Fermi level is at the mobility edge, 
Coulomb interactions between electrons in the lowest LL strongly suppress the quasiparticle 
density of states (DOS) close to the Fermi level. In certain range of energies DOS now 
exhibits a power-law dependence on the energy (at T=0) or on the temperature (at T ^ 0) 
with power 77/2, and our perturbative calculation suggests that right at the Fermi level DOS 
vanishes. This implies a weak zero-bias anomaly in the I-V characteristics of the tunneling 
into the electronic layer which would provide a direct information on the anomalous diffusion 
exponent. We discuss the experimental conditions under which this effect could be observed. 

2 Quasiparticle density of states 

Consider a 2d electron layer in a magnetic field strong enough so that all electrons are 
polarized and in the lowest LL. The LL mixing due to the disorder potential or interactions 
is neglected and we work in units in which the magnetic length I = {Tic/eB) 1 / 2 is unity. 
If ipnir) are the lowest LL eigen-f unctions of a particular realization of the disorder, the 
Hamiltonian in this basis is given by 

H = E na>t a n + ^ (m,n\v c \q,p)a+a+a p a q , (1) 

n ^ n,m,p,q 

where v c = e 2 /r is the Coulomb interaction and we measure energy from the Fermi level. 
First we assume T = 0, and calculate the self-energy correction to the single-particle Green's 
function G m (uj) = (u — E m — S m (c<j)) _1 within the perturbation theory. The real part of the 
disorder-averaged self-energy at energy E, 

E E (u) = N \E)(Y 6{E - £ m )£ m H>, (2) 

m 
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determines to the lowest order in the interaction, the quasiparticle DOS via ||10||: 



N{E) = N (E)(1 + -^ReE E (E))-\ (3) 

Nq(E) is the average DOS of non-interacting electrons in the lowest LL and in the ran- 
dom potential and angular brackets denote the averaging over disorder. Consider first the 
exchange contribution to the self-energy: 

72; 



AS, 



mL dE 'wr s(E ' E '' q)v - M (4) 

where the disorder-averaged two-particle spectral density is defined by 

S(E,E',q) = f ^oq)(^(j;^ - - i;)C(r>;(0)^(rU(0)). (5) 

n,m 

For small momentum and energies close to the mobility edge the two-particle spectral density 
takes the familiar diffusive form IHl, [111 



N (0) D(v,q)q 2 f . 

S{v, q) = 2 , n2 / — w> 6 

7r v l + D l {y, q)q^ 

and we assume the generalized diffusion constant D(i/,q) = hDf(q 2 /cNo(0)i>), v = \E — E'\, 
c and D are constants and f(y) — 1 for y < 1 and f(y) = y~ ri ^ 2 for y > 1 ||. It is essential 
for our discussion that even though we are calculating the lowest order contribution to the 
average self-energy we include the effect of screening on Coulomb interaction in eq. 4. By 
definition, the dynamically screened Coulomb interaction is given by: 

v scr (u, q) = v c (q)(l + v c (q)H(u, q))' 1 (7) 

and we assume that the polarization function has it's standard RPA form: 

N (0)D(u J ,q)q 2 

q) = — — —-, r-j. 8 

— iuj + D(u, q)q l 

As will be shown shortly, this assumption is by no means essential for our results. To 
calculate the correction to the quasiparticle DOS, we will need only the static limit of the 



screened Coulomb interaction. Thus, as long as the polarization function has a diffusive 
form, the effects beyond simple RPA approximation which would be represented by a more 
complicated diffusion constant, exactly cancel out in the calculation. Important point is that 
the static part of the screened interaction is short ranged. Using the equations 4,6,7,8 after 
some calculation we find, for E « 0: 

d p-yewm _ F exc {x,E) ,\E\ v 

dE ReEE {E) ~ 2rrg ( X } ' (9) 
where the function F exc (x, E) is given by an integral 

f 1 dt 

F exc (x,E)=X - ; r. 10 

V ' ' J(\E\/A)l/l t l -V{t + x) V ' 

9 = &xx/( e2 /h) is the dimensionless dissipative conductance a xx = e 2 DN (0), the en- 
ergy scale is A = nT/cg, F = 2hD is the half-width of the disorder-broadened LL, and 
x = (e 2 /ed)/T, where d is the average distance between electrons. The result in eqs. 9 and 

10 is the same as if we used the screened interaction in static approximation in our calcu- 
lation. Thus for the present purposes, we could have set ui = from the beginning in the 
polarization function 8 so that complications related to the form of the diffusion constant 
in that expression cancel out. Let us now analyze the equations 9 and 10. If x << 1, then 
for x 2 « \E\/A « 1 we have dReZ e £ c (E)/dE w (x/2irg(l - T]))(\E\/ A)- 1 / 2 (assuming 
r] < 1), just like the result would be if the bare Coulomb interaction was used in the calcula- 
tion. However, closer to the Fermi level, i.e., \E\/A ~ x 2 and smaller, one obtains the novel 
power-law dependence on the energy: dReT, e ^ c (E) jdE fa (l/2ngr])(\E\/ A)~ v / 2 . 

If one would perform the previous calculation using the delta-function interaction instead 
of v scr (u), q) in eq.4, close to the Fermi level self-energy would again diverge with the power 
77/2, only with a different prefactor. This means that the same term must also exist in 
the Hartree contribution to the self-energy, since the two should cancel each other if the 



interaction is infinitely short ranged [10, 11] and electrons are completely spin-polarized. 



Taking both the exchange and the direct contributions into account, we finally obtain that 



at zero temperature, sufficiently close to the Fermi level the correction to DOS reads: 

5N(E) F(x) \E\ 2 



11) 



iV (0) 2tt^ v A ' 

where F(x) = r]F exc (x,0) - F dir (x), and [1|, F dir = N (0) J d 2 qv scr (0,q)/ir. The last 



integral gives 

F dir (x) = 2x-2x 2 ln(l + aT 1 ), (12) 

and x is the ratio between the interaction and disorder energy scales as defined earlier. 

Since we found that at T = frequency dependence of the screened interaction does 
not affect the behavior of DOS sufficiently close to the Fermi level, at finite temperatures we 
may neglect it completely. The simple calculation then yields the result: 

SN(E,T) _F(x) T , \E\ 

JVo(O) "2^ ( A } t{ ^r> (13) 

for \E\,T sa 0, where 

r+oo pt+z 

™ = r itHf^+D^ (14) 



oo 



The function f(z) = z~ v ^ 2 for z >> 1 and /(0) = 0.92 (for 77 = 0.5). Thus, at zero 
temperature we recover the result in eq. 11. At T ^ 0, DOS acquires a power- law dependence 
on the temperature instead on the energy. 

A few remarks are in order at this point. The obtained power-law should be compared to 
the logarithmically divergent correction to DOS in d > 2 implied by one-parameter scaling 
relation rj = 2 — d \14\\ . The exponent rj is determined completely by the geometry of 



the extended state and is not in any sense a small parameter in the problem. The found 
power-law is therefore a truly distinct behavior from the standard logarithmic correction 
to DOS. This should be contrasted with the situations where the exponent itself is a small 
quantity so that the expansion in powers of this exponent to the leading order agrees with the 
logarithmic behavior (for a similar scenario see for instance ref. 15). Also, the assumption 
that the Fermi level is right at the mobility edge so that the diffusion constant is a function 
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of the combination q 2 /uj and not only of the momentum is important. If this was not so, 
i. e. if the Fermi level was slightly off but close to the mobility edge, the power law in 
eq. 11 would be cut off below a finite energy defined by £(0) ~ DE~^/ 2 , where £(E) is the 
localization length at energy E (E = still defines the Fermi level). For the frequencies 
smaller than the E cut diffusion constant becomes a function of combination (?£(0)), and 
independent of frequency. With 77 larger than zero the perturbative correction of the DOS 
then would show a crossover from the power law divergence as in eq. 11 for E > E cut to the 
standard logarithmic divergence for E < E cut . 

The divergence of the first-order correction to DOS in our problem indicates the break- 
down of the perturbation theory very close to the Fermi energy. Since from the eq. 10 we saw 
that the power r]/2 turns on for (\E\/A) ~ x 2 and below, one might wonder whether that 
is in the region where higher order terms are already significant. The point where the first 
order correction becomes of order unity marks the region of energies where are our lowest 
order perturbation theory becomes insufficient: this happens at (\E\/A) ~ (F(x)/2nr]g) 2 / r '. 
The function F(x) behaves like x v /(l — 77) for x < 0.01, peaks at x ~ 0.2 with the value 
of 0.26, and roughly stays constant for 0.2 < x < 1. Consider first the physically more 
relevant regime 0.2 < x < 1: then x 2 ~ 0.1 and the perturbation theory breaks down for 
E/A « 10 -3 << x 2 . Here we assumed g ~ 0.5 and 77 m 0.5 0, |5|, |6|]. The important point is 
that there is a considerable range of energies for which our simple lowest order perturbation 
theory is valid and the result is indeed given by eqs. 11 and 13 (see Figure 1). One may 
also consider the regime when x — > 0: in that case the breakdown occurs for \E\/A « x 2 , so 
that again right before it occurs the non-trivial power law appears, although the result is not 
as clear-cut as for larger x. Even though we can not say with certainty what happens very 
close to the Fermi level where the problem becomes non-perturbative, we still expect however 
that the perturbation theory is qualitatively correct and that DOS will indeed vanish right 
at the Fermi level. Standard scaling arguments would suggest that ultimately DOS still goes 
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to zero as a power-law, but with the exponent determined by the full interacting theory. 
Recent numerical calculations which include Coulomb interactions at the Hartree-Fock level 
find, quite surprisingly, both rj and the localization length exponent to be the same as in 
the non-interacting problem [lf|. This points to the intriguing possibility that the exponent 



which ultimately determines the behavior of the DOS at the Fermi level is not different from 
the one we calculated. This issue deserves more attention in the future. 

3 Tunneling 

The obtained correction of the quasiparticle DOS will modify the thermodynamic quantities 
when the Fermi level is at the mobility edge. For instance, the specific heat will acquire a low- 
temperature correction: 5C V /T oc T r, l 2 . However, thermodynamic quantities are difficult to 
measure since a typical 2d electron layer in GaAs heterostructure contains very few (~ 10 11 ) 
particles. Recently, tunneling experiments have proved to be a new useful tool for studying 
2d electron systems |17 , [T~8| , |19| . In what follows, we therefore consider a typical tunneling 



experimental setup consisting of two identical, parallel 2d electron layers separated by a 
semiconductor barrier. When a small voltage is applied between the layers, the tunneling 
current is: 

I oc N 2 (eV,T)V (15) 

with the constant of proportionality being roughly the inverse of the tunneling resistance for 
non-interacting electrons. When the Fermi levels in both layers are tuned to the mobility 
edge, the tunneling conductance should exhibit a power-law anomaly as a function of the 
voltage (T « eV) or of the temperature [eV « T) with power 77. Let us now estimate the 
relevant temperature (or voltage) scale over which the anomaly should become observable. 
DOS will start to show a power-law behavior when the correction becomes comparable to 
unity. This yields the temperature scale T pa (nT /cg)(F(x)/2%r]g) 2 / ri as already discussed. 
We assume the numbers for g pa 0.5 and rj pa 0.5 as before. The value of the constant c is 
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disorder dependent and less well known. Here, we take c ~ 60, which was found for the white- 
noise random potential ||. All the numbers put together give the relevant temperatures to 
be T oc 10 -4 * T. Note that according to our discussion, this is deep in the region where the 
power is determined by the screened Couloumb interaction. 

In a typical experiment, GaAs samples of very high mobility (for example, \x = 2 * 
10 G cm 2 V/s in ref. 18) are used. From the zero-field mobility we can estimate the LL 
broadening in the following way: assume that the scattering in the system is caused by 
the potential V(r) = A2j5(r — rj), where {fl} are 2d coordinates of uniformly distributed 
scatterers with the density n imp . The zero-field mobility is [A — eT /m, r = h/(n imp X 2 7iN ) 
is the scattering time in Born approximation and Nq = 2rmr/h 2 is 2d electron DOS in zero 
magnetic field. The band mass in GaAs is m = 0.07m e . In strong magnetic field, Born 
approximation gives the lowest LL broadening in the same potential T « (X 2 ni mp /l 2 ) l l 2 . At 
the magnetic field of 10T one then obtains V « 6K, an order of magnitude smaller than 
the Coulomb energy scale which is around 50K (at 10T, hu c ~ 160K). We treated the 
interactions as a perturbation to the disorder problem, so we need the LL broadening to be 
at least of the same order of magnitude as the Coulomb energy scale to be in the regime 
where our calculation is applicable. Presented estimate suggests that T is proportional to 
the inverse of the square root of the mobility, hence decreasing the mobility of the sample 
by a factor of 100 would make T m 60K, somewhat larger than the Coulomb energy but 
still sufficiently smaller than the cyclotron energy so that neglecting LL mixing is still a 
reasonable approximation. The temperature when the tunneling anomaly should occur is 
then ~ O.Oli^, or in terms of the voltage around lfiV. 

4 Conclusion 

To summarize, we demonstrated using the perturbation theory that the interaction effects in 
2d electron system in high magnetic field at the metal-insulator transition cause a diverging 
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power-law correction of the quasi-particle density of states close to the Fermi level. This 
is interpreted as a precursor of the power-law vanishing of the density od states at the 
Fermi energy. The power in the perturbative regime is given by the anomalous diffusion 
exponent, so a tunneling experiments may be used to test the novel diffusion behavior in 
high magnetic field discussed by Chalker. An estimate of the energy scales involved shows 
that if a double-well GaAs heterostructure of zero-field mobility of the order of 10 4 cm 2 V/s is 
used in the tunneling experiment at fields of 10T, a zero-bias anomaly in the current- voltage 
characteristics should develop for temperatures of the order of lOmK and voltages ~ 1/j,V. 
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Captions: 

Figure 1. Schematic view at the behavior of the quasiparticle density of states for small 
energies with the Fermi level at the mobility edge. The value of the small parameter x 2 pa 0.1 
marks the crossover point from the the regime where the behavior is essentially determined by 
the bare Coulomb interaction (power 1/2) into the regime where screening becomes effective 
and determines the non-trivial power law (power r]/2). In the shaded region the problem 
becomes non-perturbative. 
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